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The electron-doped silicene under the influence of the biaxial tensile strain is predicted to be the 
phonon-mediated superconductor. By using the Eliashberg formalism, we investigate the thermody¬ 
namic properties of the superconducting silicene in the case when the tension is 5% and the electron 
doping equals 3.5 x 10 14 cm' 2 . Under such conditions, silicene monolayer is expected to exhibit the 
highest superconducting transition temperature (Tc). In particular, based on the electron-phonon 
spectral function and assuming wide range of the Coulomb pseudopotential values (gP £ (0.1,0.3)) 
it is stated that the superconducting transition temperature decreases from 18.7 K to 11.6 K. Sim¬ 
ilar behavior is observed in the case of the zeroth temperature superconducting energy gap at the 
Fermi level: 2A(0) £ (6.68, 3.88) meV. Other thermodynamic parameters differ from the predictions 
of the Bardeen-Cooper-Schrieffer theory. In particular, the ratio of the energy gap to the critical 
temperature changes in the range from 4.14 to 3.87. The ratio of the specific heat jump to the 
specific heat in the normal state takes the values from 2.19 to 2.05, and the ratio of the critical 
temperature and specific heat in the normal state to the thermodynamic critical field increases from 
0.143 to 0.155. It is also determined that the maximum value of the electron effective mass equals 
2.11 of the electron band mass. 
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I. INTRODUCTION 

Recent research in the field of nanotechnology have 
led to the synthesis and characterization of various two- 
dimensional materials [Tj. The unique geometries of 
these novel structures are one of the main origins of 
their extraordinary physical and chemical properties 0, 
\2\. Among different applications, the possibility of using 
such low-dimensional systems in the domain of nanoscale 
superconducting devices is of growing interest [ 3 ], 0 , 13 - 
In this respect, the one-atom-thick two-dimensional 
form of carbon, known as graphene [6], attracted excep¬ 
tional attention in recent years, when comparing to the 
other carbon allotropes H, 0: 0- However, various 
theoretical calculations demonstrate that the phonon- 
mediated superconductivity does not occur in the in¬ 
trinsic graphene, due to the weak electron-phonon cou¬ 
pling constant m, mi Pi- This fact follows the case 
of graphite, where the induction of the superconducting 
phase is possible only via the chemical process know as 
intercalation Ba¬ 
in particular, it was suggested, by using the first- 
principle calculations, that the conventional supercon¬ 
ducting state with notable high critical temperature (Tc) 
can be observed in the hole-doped graphane (a fully 
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hydrogenated graphene) [14], |15J or in the lithium¬ 
decorated graphene H3j, [E], [H|. Due to these findings 
this direction of research can be considered as a promis¬ 
ing and important one. However, the direct experimental 
evidences are still lacking. 

Another noteworthy material for the superconducting 
nanoelectronic applications is silicene (the graphene-like 
monolayer of silicon) [18| . In general, the main advan¬ 
tage of this material is the fact that it combines some 
of the graphene intriguing properties and can be rela¬ 
tively easy incorporated into the existing silicon-based 
electronics m- Moreover, the preliminary results on the 
synthesis of silicene are already available [201 - 

From the point of view of the superconducting prop¬ 
erties, it is important that pristine silicene is character¬ 
ized by the buckled structure, which distinguish it from 
the graphene and favours stronger electron-phonon cou¬ 
pling in this material |21j . Some theoretical works, using 
random-phase-approximation (RPA), have proposed pos¬ 
sible singlet d + id' chiral supercondutivity in undoped 
bilayer silicene [22] or quantum phase transition to the 
triplet /—wave superconducting phase in doped silicene 
under a perpendicular external electric field [23] . Encour¬ 
aging are also recent experimental results which suggest 
that the induction of the superconducting state in sup¬ 
ported silicene may be possible [23] . 

On the other hand, latest theoretical investigations 
predict that the electron-doping and the influence of the 
biaxial tensile strain induce superconducting state char¬ 
acterized by the critical temperature much above 10 K 
[2Tj . In particular, for the electron doping ( n e ) equals 
3.5 x 10 14 cm -2 and tension of 5%, the analytical McMil¬ 
lan 125] formula gives Tc ~ 17 K. This outcome is promis¬ 
ing and motivates our studies. 

In the present paper, we concentrate ourselves on the 
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analysis of the superconducting phase induced in silicene 
under the conditions mentioned above. In the consid¬ 
ered case, the electron-phonon coupling constant exceeds 
the weak coupling limit (A > 0.5 [2S])- Due to this 
fact we conduct our calculations within the framework 
of the Eliashberg formalism [27], which allows us to cal¬ 
culate the thermodynamic properties of the supercon¬ 
ducting phase at the quantitative level. Our calculations 
based on the electron-phonon spectral function a 2 F(uj) 
obtained in [21j . by using the density functional theory 
within the local-density approximation. The numerical 
methods adopted in the present work based on the self- 
consistent iterative procedure presented in m,m, wi¬ 


ll. THE NUMERICAL AND ANALYTICAL 
RESULTS 

In order to compute all thermodynamic properties of 
interest, the Eliashberg equations are solved on the imag¬ 
inary axis and in the mixed representation (defined simul¬ 
taneously on the imaginary and real axis). The stability 
of the solutions in both cases is achieved in the tem¬ 
perature range from To = 1 K to Tc, for the assumed 
1100 Matsubara frequencies: w m = ^(2m — 1), where 
/3 = 1/fceT, and fcs denotes the Boltzmann constant. 

The Coulomb pseudopotential (/x*) models the de¬ 
pairing interaction between the electrons and beside the 
Eliashberg function is the second input parameter in the 
Eliashberg equations. In fact there are two well-known 
experimental methods to determine fi*. One is based on 
the isotope effect [25], second is based on the inversion of 
tunnelling data J3D|. It should be emphasized, that the 
isotope effect can be used only when a set of isotopes is 
available and the tunnelling experiments require strong 
or medium coupling superconductors in order to give suf¬ 
ficient structure in the current-voltage characteristic ED- 
The physical value of Coulomb pseudopotential is hard 
to calculate from first-principles, so it is usually treated 
as a material-dependent adjustable parameter chosen, for 
examples within the framework of the Eliashberg formal¬ 
ism, to reproduce the experimental value of critical tem¬ 
perature [32] . However, due to absence of the experimen¬ 
tal results for silicene, a wide range of the Coulomb pseu¬ 
dopotential values is taken into account: jj* G (0.1, 0.3). 

In figure [l] (A), the superconducting order parameter 
on the imaginary axis for selected values of the tempera¬ 
ture and the Coulomb pseudopotential is presented. The 
maximum value of the order parameter (A m= i) decreases 
with the growth of T and /j,* . This fact can be clearly ob¬ 
served in figure[T](B) where A m= i (T) function is shown. 
On the basis of these results, we note that the A m=1 val¬ 
ues can be well parameterized by the following formula: 


A m=1 (T, M *) = A m=1 (/x*)yi-^ , (1) 

where: A m= i (/i*) = 18.30 (/x*) 2 — 14.08/i*+ 4.51 and the 


fitting parameter a is equal to 3.4. 




FIG. 1: (A) The dependence of the order parameter on the 
number m for the selected temperatures and Coulomb pseu¬ 
dopotential values. (B) The full dependence of the maximum 
value of the order parameter on the temperature for selected 
/+. (C) The critical temperature as a function of the Coulomb 
pseudopotential. The circles correspond to the exact numer¬ 
ical solutions of the Eliashberg equations. The triangles and 
squares represent the results obtained using the Alien-Dynes 
and McMillan formula, respectively. 


The superconducting transition temperature is defined 
as the temperature at which the order parameter van¬ 
ishes: A m= i (Tc,/x*) = 0. In the case of silicene it is 
stated that Tc is high in the entire range of the Coulomb 
pseudopotential values. In particular, Tc G (18.7,11.6) 
K for /x* G (0.1,0.3). It should be noted that the max¬ 
imum value of the critical temperature for /x* = 0.1 de¬ 
termined by us is significantly higher then the value pre¬ 
dicted in paper [2D, where [Tc]™=qa = 16.40 K. This 
situation is caused by the fact that in paper [2D the su¬ 
perconducting transition temperature was estimated by 
using McMillan formula [25]: 


, rp w ln 

knTc = — exp 


-1.04(1 +A) 

A - n* (1 + 0.62A) 


( 2 ) 
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where uq n is the logarithmic averaged phonon frequency 
and equals 18.52 meV for the tension of 5%. 

A full dependence of Tc on p* is plotted in figure [l] 
(C). Presented results are obtained using the Eliashberg 
formalism, Alien-Dynes expression m and the McMil¬ 
lan formula [25y It can be observed that the analyti¬ 
cal results underestimate the critical temperature, espe¬ 
cially for the high values of the Coulomb pseudopoten¬ 
tial. Moreover, the Alien-Dynes expression much better 
predicts Tc than the McMillan formula. 

In figure [2] we present the results for the wave func¬ 
tion renormalization factor. The identical values of tem¬ 
perature and Coulomb pseudopotential as for the order 
parameter are chosen. It is found that the function 
Z m= i(T) increases together with the temperature and 
the Coulomb pseudopotential value. 



FIG. 2: (left panel) The form of the wave function renor¬ 
malization factor on the imaginary axis for selected values of 
temperature and Coulomb pseudopotential, and (right panel) 
the dependence of the maximum value of the wave function 
renormalization factor on the temperature. 

The values of the function Z m —i (T,p*) can be esti¬ 
mated by the formula: 

Zm=i(T,n*) = Z m —i(fj, *) (3) 

+ [Z m -1 ( T c ) - Z m= i (p*)] ^ —^ , 

where: Z m —i (p *) = —0.48 (p*) 2 + 0.35 p* + 1.98 and 
Z m —i (Tc) = 1 + A. Note that for T = Tc, the max¬ 
imum value of the wave function renormalization factor 
is independent of p*. 

The temperature dependence of A m= i is reflected in 
measurable thermodynamic parameters as the specific 
heat C S (T) of the superconducting state or the ther¬ 
modynamic critical field Hc(T). In order to investigate 


these properties from the solution of Eliashberg equa¬ 
tions on imaginary axis, we have evaluated numerically 
the free energy difference between the superconducting 
and the normal state (A F = F s — F N ). This is given by 
the formula [33]: 


A F 

JW) 


x 


M 

-^£(7^TaI' 

“ n= 1 




-z. 


N 


\LO n 


V U n + A n 


(4) 


where p ( 0 ) is the value of the electron density of states 
at the Fermi energy level. Symbols Z% and Z^ de¬ 
note the wave function renormalization factors for the 
superconducting (A m= i ^ 0) and for the normal state 
(A m= i = 0), respectively. 

The thermodynamic critical field and deviation func¬ 
tion of the thermodynamic critical field are calculated 
from the free energy difference: 



(5) 

T) \l ( T ) 2 ] 

H c ( 0) L \T C J \ 

(6) 

The specific heat difference AC = C s — C N 
to A F through: 

is related 

AC (T) 1 d 2 [A F/p (0)] 

W( 0) “ P d (k B T) 2 ’ 

(7) 


while the specific heat in the normal state is defined as: 
C N = 7 //J, where 7 is the Sommerfeld constant: 7 = 
(2/3)7t 2 (1 + A) UbP (0). 

The results obtained for silicene under the tension of 
5% are presented in the lower panel of figure [3] (A). From 
the physical point of view, the negative values of A F/p{ 0 ) 
inform us about the thermodynamic stability of the su¬ 
perconducting state in the temperature range from To 
to Tc- It can be seen that the increase of the Coulomb 
pseudopotential substantially weakens the stability of the 
superconducting state in silicene. 

The upper panel in figure [3] (A) presents the influ¬ 
ence of the temperature on the ratio He / \Jp (0). It can 
be seen that the thermodynamic critical field decreases 
with the increasing temperature, taking the zero value for 
T = Tc- Let us notice that the maximum values of the 
considered function equal: He (0) /\Jp (0) € (15.8,9.4) 
meV for p* £ (0.1,0.3), where He (0) = He (To). 

In figure [3] (B), the specific heat for the normal and 
superconducting state is presented. We can note that 
with the increase of the temperature the specific heat 
of the superconducting state grows strongly, reaching its 
maximum at the critical temperature. The characteristic 
specific heat jump at the critical temperature is marked 
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with the vertical line. We can observe that the value of 
the jump , just like the value of the thermodynamic criti¬ 
cal field, decreases with the growth of the Coulomb pseu¬ 
dopotential. In particular this fact can be described by 
the following ratio: [AC (T c )] * =0 3 / [AC (T C )] M , =0 1 = 
0.58. 

In figure [3] (C) we supplement our results with the cal¬ 
culated thermodynamic critical field deviation as a func¬ 
tion of the temperature for selected values of Coulomb 
pseudopotential. The positive values of D function cor¬ 
respond to the strong electron-phonon coupling (A > 1) 
and D is negative for the weak coupling limit (A < 1) [53] , 
Our results confirm that for the investigated material the 
electron-phonon coupling is strong. Results predicted by 
the Bardeen-Cooper-Schrieffer theory (BCS) [35], [36] are 
presented with a dashed line Ell¬ 
in the next step, on the basis of the calculated 
thermodynamic functions, we determined the dimen¬ 
sionless ratios Rc = AC (Tc) /C N (Tc) and Rh = 



T/T c 


FIG. 3: (A) The free energy difference (lower panel) and the 
thermodynamic critical field (upper panel) as a function of 
temperature. (B) The specific heat of the superconducting 
and normal state as a function of temperature for selected 
values of Coulomb pseudopotential. (C) Critical field devia¬ 
tion as a function of the temperature. 


T c C n (Tq) /Hq (0). In the framework of the BCS model 
these parameters have the universal values: [-RcJbcs = 
1.43 and [i?#] BCS = 0.168 [35], [36]. Taking into account 
obtained results, we found that the values of Rc and Rh 
for silicene under the tension of 5% significantly diverge 
from the values predicted by the BCS model, in particu¬ 
lar: R c G (2.19,2.05) and R H G (0.143,0.155). 

The Eliashberg equations in the mixed representation 
are next solved for the identical range of the tempera¬ 
tures and the Coulomb pseudopotential as the Eliashberg 
equations on the imaginary axis. 

On the basis of the solutions of the Eliashberg equa¬ 
tions in the mixed representation the exact value of the 
order parameter can be obtained from the relation [38 .: 

A (T) = Re [A (w = A (T ), Tj] . (8) 

The zeroth temperature superconducting energy gap 
at the Fermi level for the investigated case of silicene 
takes the following values: 2A(0) = 2A(T 0 ) G (6.68,3.88) 
meV for fj* G (0.1,0.3). The most interesting from 
the physical point of view is dimensionless ratio R/\ = 
2A(0)/ksTc- In the framework of the BCS theory: 
[I?a]bcs = 3-53 [35], whereas in our case, even for a large 
values of /i*, the ratio significantly exceeds the values pre¬ 
dicted by the BCS model, in particular R& G (4.14,3.87). 
This situation is caused by the fact that in the BCS 
model, the strong-coupling and retardation effects are 
not taken into account. Therefore this theory is valid only 
for values of the gap functions which are small compared 
to the Debye frequency. Strong-coupling and retardation 
effects are taken into consideration in Eliashberg theory 
and can be characterized by the ratio k B Tc /^i n - In the 
weak-coupling limit, one can assume: [k B Tc /uj\ n ] BCS — I 
0. For two-dimensional silicene layer under the ten¬ 
sion of 5% we obtained: [k B Tc/oj i n ] M * = o 1 — 0.087 and 
[k B T c /iv i„] m * =0 . 3 ^ 0.054. 

Moreover, the order parameter function on the real 
axis allows to calculate the total normalized density of 
states [271 : 


NDOS (w) 


DOS s (w) _ 
DOS n M “ 


|w — *r| 

y/(u,-tT) 2 -A* M 


(9) 

where the symbols DOSs ( w ) and DOSn (w) denote the 
density in the superconducting and normal state, respec¬ 
tively. The pair breaking parameter T is equal to 0.15 
meV. 

The calculated total normalized density of states 
for different temperatures, and for selected values of 
Coulomb pseudopotential is presented in figure [4] These 
results allows us to trace the size of the superconduct¬ 
ing gap with increasing values of the temperature and 
Coulomb pseudopotential. The symmetric maximas can 
be clearly observed for T < Tc- It is worth empha¬ 
sizing, that the superconducting gap between the two 
symmetric maximas, is much larger than in the lithium¬ 
decorated graphene (LiCg) and is comparable with L^Cg 
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FIG. 4: The dependence of the total normalized density of states on the frequency for selected temperature and values of the 
Coulomb pseudopotential. 


[25] , Above the critical temperature, the gap is no longer 
visible and silicene can reveals metallic behavior. 

In the last step, the wave function renormalization fac¬ 
tor on the real axis (Z (cu)) is determined. In the frame¬ 
work of the Eliashberg formalism the real part of Z ( oj ) 
enables the determination of the electron effective mass 
(to*). In particular, the ratio of to* to the electron band 
mass ( m e ) is given by: to*/ m e = Re [Z (0)]. Let us no¬ 
tice that R e[Z (0)], similarly as Z m —\, takes the highest 
value for T = Tq. Thus, [m/] max is equal to 2.11m e . 

III. SUMMARY 

Using the Eliashberg approach, we have studied the 
behaviour of the superconducting critical temperature, 
energy gap, free energy difference between the supercon¬ 
ducting and normal state, thermodynamic critical field 
and the specific heat in a strongly coupled electron- 
doped silicene under the tension of 5%. The Coulomb 
pseudopotential effects on the thermodynamic properties 
have been extensively studied based on electron-phonon 


spectral function. For /z* £ (0.1,0.3) the critical tem¬ 
perature and zeroth temperature superconducting energy 
gap at the Fermi level decrease from 18.7 K to 11.6 K and 
from 6.68 meV to 3.88 meV, respectively. 

Other thermodynamic parameters differ from the pre¬ 
dictions of the BCS theory. In particular, for the dimen¬ 
sionless ratios of the calculated thermodynamic functions 
we reported the following results: Ra £ (4.14,3.87), 
Rc £ (2.19,2.05) and Rh £ (0.143,0.155). It is con¬ 
nected with fact that the Eliashberg formalism in con¬ 
trast to BCS model does not omit the strong-coupling 
and retardation effects. 
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